The effect of entropy on macroions adsorption 
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We study macroion adsorption on planar surfaces, through a simple model. The importance of 
entropy in the interfacial phenomena is stressed. Our results are in qualitative agreement with 
available computer simulations and experimental results on charge reversal and self-assembling at 
interfaces. 
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The physics of macroions adsorption at charged inter- 
faces is subject of recent theoretical and experimental 
studies 0, S S IE @ because its potential technologi- 
cal applications. It is known that in the neighborhood of 
a charged surface, low-concentrated solutions of monova- 
lent ions produce an exponentially decaying charge dis- 
tribution which is well described by Poisson-Boltzmann 
equation. However, adsorption of multivalent ions dis- 
plays important deviations from the Poisson-Boltzmann 
picture and there is not yet a fully accepted theory which 
well accounts for this situation 

Bii 

In inhomogcneous 

charged fluids, charge reversal (CR) and charge inver- 
sion (CI) have been reported as mechanisms with which 
systems minimize free energy. CR is the overcompen- 
sation of the surface's charge by a contiguous layer of 
oppositely charged ions from the fluid. Because of elec- 
trostatic equilibrium, next to the CR layer, a second layer 
of ions (with same charge sign as that of the surface) is 
formed, producing a CI In molecular engineer- 

ing, CR and CI are basic mechanisms for self-assembling 
polyelectrolyte layers on a charged substrate 0, reverse 
mobility experiments , self-assembled DNA-lipid mem- 
brane complexes |lfj L anomalous macroions adsorption 
on a Lagmuir film [2j, and novel colloids stabilization 
mechanisms 0. 

Although CR was reported since the early 80's [Tllll^] 
and some advances in its understanding have been done, 
the mechanism has not been completely understood. 
Some authors recognized the relevance of electrostatic 
and ionic size correlations in CR 0, @, 0] . In particu- 
lar, the unsymmetrical size effects (between small coions 
and counterions) has been considered by Greberg and 
Kjellander 0. Their analysis of CR and CI is based 
on mechanistic arguments and is basically correct, al- 
though incomplete: the role of the ionic concentration, 
entropy and charge asymmetry is not discussed, which 
are not minor effects 6], as we will see. Recently CR by 
macroions of an oppositely charged surface was consid- 
ered by Shklovskii et al. [l3j- They assert in proposing 
that, when macroions overscreen a charged surface, these 
macroions are highly ordered in a bidimensional Wigner 
crystal (WC) structure, as a configuration of minimal 



free energy. However, because in the WC theory ionic 
correlations are only considered at the level of a two 
dimensional electrons liquid at zero temperature, this 
approach leads to incorrect conclusions. For example, 
it is predicted that CR decreases by increasing surface 
charge density (see Fig. 3 in ref. |13|) whereas molecu- 
lar simulations show the opposite |lJ,lli|. Nevertheless, 
WC theory has a good agreement with non-zero Monte 
Carlo simulations of a two-dimensional electron liquid 
(lif and zero temperature molecular dynamics calcula- 
tions of a two-dimensional layer of adsorbed ions [l7| . 
The role of entropy in the attraction of like-charged par- 
ticles has been recognized in the past Here we will 
analyze the effect of electrostatic and short range correla- 
tions on macroions adsorption at interfaces. A decreasing 
of the system entropy always produces an increasing of 
macroions adsorption. Thus, it is concluded that self- 
assembling at interfaces, CR, CI and like-charged attrac- 
tion are mainly ruled by entropy. 

For the present study we consider a uniformly charged 
planar surface with oo charges per unit area contigu- 
ous to a macroions solution. In our model, we in- 
clude the two relevant forces of the system: the long 
range coulombic interaction and the hard-core interac- 
tion forces. Therefore macroions are considered to be 
charged hard spheres (with charge Q = Ze, diameter 
a and number density p\ being e the electronic charge) 
plus an additional amount (p c ) of neutralizing counteri- 
ons (with charge q — z c e and radius rg — a c /2). The 
mixture is embedded in a continuous media of dielectric 
constant e, set equal to that of the charged surface in 
order to avoid image charges. By the electroneutrality 
condition -a = | °° p e i(x)dx with p ei = £) i=1 (sc) , 
and pi (x) being the local concentration profile of the i-th 
species. Our model is an extension to that considered 
by Kjellander and Greberg @. The WC theory [3] has 
some similarity with our model. However, whereas we 
use direct (macroion-macroion, macroion-counterion and 
macroion- wall) coulombic interactions, in WC screened 
interactions are used. More importantly, short range cor- 
relations and entropy (of fundamental importance for fi- 
nite temperature calculations) are included in our theory 



Typeset by REVTpjX 



2 




20 40 60 80 100 120 20 40 60 80 100 120 

x[rJ x[r ] 



FIG. 1: Reduced density profile g(x) for negative macroions 
(Z — 20) next to an positively charged surface (a a = 
0.03C/m 2 ). We present results for n = 0.002, 0.01,0.05. In 
the inset we show the corresponding adsorbed charge densities 
a(x). 




We study the system by means of integral equations 
using the hybrid closure hypernetted chain/mean spheri- 
cal approximation (HNC/MSA). This approximation, for 
simpler models of charged inhomogeneous fluids, has 
shown to be in agreement with computer simulations 
where CR and other effects are studied [13. Ilif . For the 
present model, we found qualitative agreement with in- 
tegral equations results and molecular simulations Q . 
The HNC/MSA equations can be derived in an straight- 
forward manner |2fJ and are given by 

gi(x) = exp |-/3u Qi (x) + ^Pj J h j (y)c ij (s)dv 3 ^j , 

(1) 

where pi is the bulk concentration of the i-th species, 
hj(x) = Qj(x) — 1 and Cjj(s) are the total and direct cor- 
relation functions, respectively; u a i(x) is the ion-surface 
interaction potential, (3 = l/k^T with fc^ the Boltz- 
mann constant and T the absolute temperature. From 
HNC/MSA, the concentration profiles (pi(x) = pigi(x)) 
are obtained directly. These are related to the ion-surface 
potential of mean force, u>i(x), by gi(x) = cxp{— f3wi(x)}. 
Eq. n can be written as Wi(x) = —ksThi[gi(x)] — 
tjji(x) + Ji(x), where ipi( x ) is the mean electrostatic po- 
tential and Ji{x) are the short-range contributions. Since 
both, ip(x) and Ji(x) are functionals of p e i(x) and are in a 
non-linear equation, the charge and size correlations are 
non-independent. In our theory macroions adsorption is 
related to both: charge and ionic size correlations. In 
the Poisson-Boltzmann description the short range cor- 
relations are not considered in any way (i. e., Ji{x) = 0). 
Therefore, PB does not describe correctly many of the 



FIG. 2: Same as Fig.[T]but a = -0.03C/m 2 . 



adsorption phenomena. 

In a system at constant volume V, fixed par- 
ticle number N and T, its entropy is given by 
S = S id + k(^§§^) NV , where 5 id is the ideal 
gas entropy and the configurational integral Zjq = 
J ... J exp {— [3Un} dri, dr^, being Un the potential 
energy of the N interacting particles. The exact compu- 
tation of S is not an easy task. However, from the expres- 
sions for s and Zm is clear that: (i) the accessible volume 
|2~l| decreases as the volume fraction (t]t = ? J2i Pi a i) 
increases. Thus, the larger tjt the lower the entropy, (ii) 
Electrostatic interactions between like-charged particles 
also reduce the system number of accessible configura- 
tions (hence entropy) by increasing the strength of the 
coulombic coupling interaction £ = /3qf/ea,i. Less en- 
tropy always implies higher order, which in an inhomo- 
geneous fluid translates in higher adsorption, as we will 
show. In this picture it is plausible charge reversal by 
monovalent ions, macroions adsorption on uncharged (or 
like-charged) surfaces, as well as the strong dependence 
of CR on r\. These effects can not be explained only with 
energy arguments. 

We will discuss our results in terms of gi(x), the local 

X 

net charge per unit area a(x) = a + J p e i(y)dy and the 



absolute maximum of this function (a* = max{cr(x) : 

x € [0,oo)}); with the charge profile given by p e i(x) = 
2 

^2 ez m p m gm(x). The effective electrostatic force on an 

TCI — 1 

ion is ff(x) = =^-qi<j{x). Therefore, cr{x) not only quanti- 
fies charge adsorption but also the wall-particle effective 
electrical force. We have solved HNC/MSA for several 
values of Z and p. In all cases the little ions (counte- 
rions) are monovalent, i.e., \z c \ = 1. In all our calcula- 
tions we have used fixed values of a = 4nm, a c — 0.4nm, 
T = 298K and e = 78.5. We present results for macroions 
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next to positively charged, negatively charged and un- 
charged surfaces. We study the influence on macroions 
adsorption of <To, Z, 77. 

In Fig. ^ we show the macroions reduced concentra- 
tion profile (RCP) for negative macroions (Z = 20, i.e., 
Q = 20e < 0) next to a positively charged surface 
((To = 0.03C/m 2 ) for three different macroions concentra- 
tions (77 = 0.002,0.01,0.05). The macroion RCP shows 
a very strong adsorption to the wall and a considerable 
amount of counterions surrounding macroions: In all fig- 
ures, the counterions RCPs are not plotted for clarity. 
A second layer of macroions is adsorbed mediated by 
counterions. In the inset we show the corresponding 
o~(x) functions: At a distance of one macroion radius 
a deep minimum is observed, corresponding to a very 
strong CR. For rj = 0.05, the maximum located around 
x = 35ro shows a CI. In the inset, the maximum located 
at x = a/2 indicates that small positive ions are adsorbed 
to the wall. Hence, positive charge is accumulated next 
to the positively charged surface. We denominate this 
effect as overcharging. This effect is energetically unfa- 
vorable, thus, it must be attributed to the short range 
correlations since it increases with 77. This prediction 
agrees qualitatively with computer simulation of Tanaka 
and Grosberg [jj. The effective wall electrical field, which 
is proportional to a(x), is positive before the first layer 
of macroions and then negative, before the second layer. 
Therefore, the electrical force is first attractive and then 
repulsive to negative ions. The behavior of the total force 
on an ion of species i, however, is implicit in the RCP, i.e., 
a gi(x) above (below) 1 implies that Fi(x) is attractive 
(repulsive) . 

Some of the physics of macroions adsorption is not ex- 
plicit when considering macroions adsorbed on oppositely 
charged surfaces. In Fig. [2 we show the macroions RCP 
for negative macroions (Z = 20) next to a negatively 
charged surface (ctq = — 0.03C/m 2 ). Counterions adsorb 
on the charged surface (not shown) and, for sufficiently 
high values of ryr, a layer of macroions is adsorbed to the 
wall, after a relatively large void of macroions. Hence 
a surface CR is induced by little ions (according to the 
o~(x) functions in the inset) and they, in turn, induce 
macroions adsorption and the formation of the macroions 
layer. For sufficiently low values of 77 = f-pa 3 there is 
no macroions adsorption. The qualitative agreement of 
our Fig. 12 results for 77 = 0.01, with Monte Carlo data 
reported in ref. Q|, is encouraging. From Figs. ^ an d 
it is seen that higher t?t produces higher adsorption. 
This prediction is in agreement with computer simula- 
tions HQ Ell- 
in Fig. |3 the RCP are shown for negative macroions 
next to an uncharged surface. For this macroions volume 
fraction (77 = 0.05) we observe colloids adsorption even 
for low charged colloids [as low as Z = 1 (not shown)]. 
Since the wall is uncharged, the strong macroions ad- 
sorption and wall overcharging is due to the low system 




FIG. 3: Reduced density profile g(x) for macroions next to 
a uncharged surface. We present results for Z = 20, 30, 50; 
77 = 0.05. In the inset are the corresponding adsorbed charge 
densities a(x). 



entropy, which in this case is due to both: high 77 and 
high £. Increasing Z increases (£), and thus decreases 
the entropy and increases macroions adsorption and in- 
ducing a surface charge. This induced charge is in agree- 
ment with the nanoparticles halo reported by Tohver [|J , 
which seems to be a new colloidal stabilization mecha- 
nism. Our results show that there is not a minimum sur- 
face or macroions charge for the occurrence of CR. There- 
fore, the assembled nanostructures at interfaces are sta- 
bilized by entropy-induced attractive forces, rather than 
the surface-macroions electrostatic correlations as it is 
proposed in the the WC theory. These findings are in 
agreement with experiments on polyclectrolyte adsorp- 
tion where it is found that CR of a charged substrate "is 
more a property of the polymer than a property of the 
surface" 0, and with previous theoretical results, where 
the importance of particle's excluded volume (entropy) 
is recognized as an attractive force mechanism between 
like charged particles [l^. Please notice that CR in Fig. [21 
through Fig- 01 is of the same order of magnitude whether 
the substrate charge is positive, negative or zero. 

In Fig. 01 a*/o-Q is plotted as a function of Z for 



r] = 0.002,0.01,0.05, for both a 



-0.03C/m 2 and 



cr = 0.03C/m 2 (in the inset). For a = -0.03C/m 2 , 
macroions and surface are both negatively charged, and 
CR is carried out by the small counterions (see discus- 
sion of Fig. |2J). Nevertheless a* /(To strongly depends 
on rj, being negligible for low rj. For o~o < 0, <7*/<to 
monotonically increases as Z with a positive concavity 
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\^~ a W > °) whereas for a Q > — w 
sults imply that oppositely charged macroions are much 
more efficiently adsorbed into a charged surface. Our 
results for oppositely-charged macroions resemble those 
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< 0. These re- 
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higher DNA adsorption on the membranes, i.e., CR on 
the membrane. This effect can not be explained by en- 
ergy arguments, but according to our findings an increase 
in macroion concentration produces higher adsorption. 

In summary, a model for macroions solution next a 
charged surface has been considered. We found that en- 
tropy plays a fundamental role in macroions adsorption 
phenomena and, in consequence, for assembly of orga- 
nized nanostructures at interfaces. Some of our predic- 
tions seem to be in agreement with previous theoretical 
studies of polyelectrolyte solutions [18j, computer simu- 
lations 0,0)^3 an d experimental results [H13. Irl Ho| . 

We gratefully acknowledge the financial support of IN- 
DUSTRIES NEGROMEX. 



FIG. 4: a* /o~o as a function of the, negative, macroions va- 
lence Z with cto = -0.03 , for 77 = 0.002, 0.01, 0.05, 0.4. In the 
inset we show the results when cto = 0.03C/m 2 . 

shown in Fig. 3 of ref. [H (WC theory) where a*/a Q 
increases in the same way with Z and ^^52%- > 0. How- 
ever, the role of entropy is not properly considered in the 
WC theory m WC there is a minimum in 00 and 
Z must be large to have overcharging which is in opposi- 
tion to computer simulations [l3 | and our theory. When 
entropy is considered we observe the following: i) a* / ' <jq 
depends more strongly on 77 than on Z 22] . ii) For a 
sufficiently high value of r\ our theory predicts a* > 0, 
even for Z — 1 , independently of the value of oq (even for 
(7 = 0, see Fig. 0}. iii) Although not shown explicitly, 
we found that a* /<jq (i.e., CR) increases by increasing 
|<7o| (positive or negative, not shown), which is in agree- 
ment with recent simulations results \l4 . 1 1 5| and, as we 
pointed out above, in opposition to the WC theory. 

From Figs. El to 0] it is clear that macroions adsorp- 
tion strongly depends on r]T- This effect is far more im- 
portant for CR than surface or ionic charge. In an ex- 
periment of macroion adsorption on a langmuir film of 
amphiphilic molecules Rondelez Q finds that macroions 
adsorption is increased (hence CR increases) as the lipid 
layer charge density decreases, by increasing the sur- 
face area. This is consistent with our results since an 
increase of the area per molecule decreases the surface 
charge density but not the macroions adsorption which 
is strongly driven by macroions excluded volume. Ac- 
cording to this the CR should also increase by increasing 
the macroions size and/or concentration. In another im- 
portant experiment Radler et al. find that the effective 
charge on DNA-cationic liposome complexes depends on 
the DNA-to-cationic lipid mixing ratio, i.e., for a cer- 
tain DNA-to-DNA distance, dDNA=di SO the complex is 
neutral. However, if the DNA solution concentration is 
increased (Idna becomes lower than di SO , which implies 
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